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STABILITY OF AN ELECTRON BEAM INJECTED INTO SPACE 


By R. Krishan Jaggi 
Manned Spacecraft Center 

SUMMARY 


This paper is a study of the stability of an electron beam injected into space, 
with the beam assumed to be confined to a cylindrical region and assumed to be in 
equilibrium under its electric charge and its thermal pressure. A dispersion equa- 
tion, derived for surface perturbations symmetrical about the axis of the beam, shows 
that the beam is stable. 


INTRODUCTION 


The subject of this paper is the stability of a beam of energetic electrons in- 
jected into the ionosphere. A beam that consisted of 10-keV electrons was fired from 
a rocket at an altitude of approximately 270 kilometers above Wallops Island, Virginia 
(ref. 1). The beam was pulsed from an accelerator to yield pulses of 0. 1- and 
1- second durations. The auroras produced by the 1- second pulses were photographed 
by ground-based cameras. From the length of the observed rays and from the fact 
that the beam power penetrated the atmosphere to an altitude of approximately 
100 kilometers, it may be concluded that there were no serious instabilities in the 
beam. 


Some of the instabilities of an electron beam have been discussed by Tidman 
(ref. 2), Thompson (ref. 3), and Jones (ref. 4). Tidman studied the instability of 
whistler modes by using a Maxwellian distribution for ambient electrons with a hump 
caused by the high-energy beam. Thompson studied the instabilities in electrostatic 
modes propagating at an angle to the magnetic field. In all the aforementioned work, 
the beam charge was not taken into account. Furthermore, the beam was assumed to 
occupy infinite space. 

Although the previously mentioned instabilities of an electron beam are impor- 
tant, it would be premature to conclude that the beam will be destroyed, because the 
electric field produced by the beam charge may suppress some of the instabilities. 

In this paper, the beam is considered to be confined to a cylindrical region of radius 
Tq. The problem to be studied is the stability of perturbations on the surface of such 

a beam. The electrostatic and whistler wave instabilities that may arise from the 
interaction of the beam with the ionosphere can produce a spectrum of perturbations 
in the beam surface. 
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SYMBOLS 

constants of integration 

magnetic field vector 

velocity of light 

electric field 

electric field vector 

electric field at the beam surface 

r and z components of E inside the beam 

r and z components of E outside the beam 

electron charge 

zero-order Bessel function of the first-kind, imaginary 
argument 

first-order Bessel function of the first-kind, imaginary 
argument 


zero-order Bessel function of the second-kind, imaginary 
argument 

first-order Bessel function of the second-kind, imaginary 
argument 

wave number 

azimuthal wave number 
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m 

electron mass 

n 

electron density of the beam 

n’ 

electron density outside the beam 

- n + n T 
n " 2 


p o 

pressure of the electron gas at the beam axis 

P 1 

pressure perturbation inside the beam 

P l’ 

pressure perturbation outside the beam 

P e 

pressure of the electron gas 

r 

perturbed radius 

r o 

unperturbed radius 

t 

time 

U 

macroscopic velocity of the electron gas 

u , u 
r z 

r and z components of u inside the beam 

u \ u ’ 
r ' z 

r and z components of u outside the beam 

z 

coordinate along the beam axis 

5 

amplitude of the perturbation 

O 

azimuthal coordinate 

l = yfl 


<Z> 

electric potential 

03 

perturbation frequency 




= ^4n 


4irne^/ m 
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on the two sides of the beam 


£Dp' = ^4nn'e^/ m 

[x] the difference between x 

surface 


Operator: 
div( ) 
v( ) 

Curl( ) 


divergence 
gradient 
V X ( ) 

AN ASSUMED EQUILIBRIUM 


It is assumed that the cylindrical beam is in equilibrium under its electric field 
and its thermal pressure. If the electron density n is uniform inside radius r_, 
the electric field is radial and is given by 


E = -E 


0 — for r =s r Q 


( 1 ) 


and 


E = -E 


0 -forrar 0 


( 2 ) 


where E^ = 2irner . For equilibrium, if it is assumed that the density of the ambient 
medium inside r = r^ is negligible, the pressure p g should vary as 


P 


e 




( 3 ) 


2 2 2 

where p^ = nn e r^ . Outside the beam, it is assumed that the ambient medium is 

at zero temperature. These assumptions are made for mathematical simplicity, and 
the stability of the configuration is not expected to change if the outside medium is 
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assumed to be at a definite temperature and if a positive surface charge is 
at r» in order to satisfy the pressure balance. In the actual case, the beamd 


served to expand. Assumption of a cylindrical beam in equilibrium at radius r^ 

implies that the expansion velocity is small compared to the phase velocity of the 
perturbations. 


THE DISPERSION EQUATION 


The equations governing the motion of electrons are 


du 

nm -jr = -neE - vP 
dt e 

( 4 ) 

div u = 0 

( 5 ) 



div E = -47me 

(6) 

-»■ 47T 0U 


Curl Curl E = — g ne 

( 7 ) 


c 


Use of equation (5) assumes the electron gas is incompressible. This assumption is 
justified in the present case if the phase velocity of the perturbations is small com- 
pared to the speed of sound in the beam. On physical grounds, the compressibility is 
expected to be important for radial pulsations that cannot take place in an incompress- 
ible medium. For the type of perturbations studied, compressibility is not required. 
Equation (7) is derived from Curl B = -(477/c)neu and Curl E = - (l/c)(0B/0t). In 
equation (7), the displacement current is neglected, which imposes the restriction of 
using only waves with phase velocity less than the velocity of light. The perturbations 
in the beam surface may be represented by 


K Z m t+Zkz+ZlO 
r = r„ + oe 


( 8 ) 


The present analysis is confined to the symmetrical case of 1 = 0. 
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If the r and z components of the perturbation velocity and electric field are 

denoted by u , u , E , and E (fig. 1) and if the perturbation in the pressure is 
r z r z 

denoted by Pj, the solution of the linearized equations can be written as 


Pj = CI 0 (kr) 


r<r 0 

(9a) 

u z = AVKr) 


r<r o 

(9b) 

u r = X A V Kr ) 


r < r o 

(9c) 

E z = ' ne kCI 0^ r ^ + nm mAI 0 (Kr) 


(9d) 

E = - ~ CI..(kr) - 
r ne 1 

n ”^AI 1 (Kr)l 


(9e) 


where K = >r 2 + v° P /^ 2 ) ; ® P = «J 4nne/m; Iq and 1^ are zero- and first-order 

Bessel functions of the imaginary argument, respectively; and C and A are con- 
stants of integration. If, for the outside medium, n' denotes the electron density, 
then the solution of the linearized equations can be written as 


Pi’ 

= GK Q (kr) 

(10a) 

u T 
z 

= DK Q (K’r) 

(10b) 

u r' 

= ^ K x (K’r) 

(10c) 

E ' 
z 

= - ^[kGK 0 (kr) + n'm m DK 0 (K’r) 

(10d) 

E r’ 

- Hl[ GK l (kr) - DK 1<H 

(10e) 
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Ambient quantities 

Perturbed quantities 

Outside 

Density - n 



Pressure = 

Pressure = 0 



Electric field = , 0, E z ‘^ 

^radial 



Velocity 0, u z '^ 

Inside 

Density - n 



Pressure = p^ 

Pressure — Pq j^l - y / r ^ 

)] 


Electric field = ^E f/ 0, 

m 

& 

Q_ 

CD* 

II 

m 

o 

N 



Velocity = ^u,, 0, u^ 



Figure 1. - Ambient and perturbed quantities inside and outside 

the beam plasma. 


where K' = 


k 2 + 


M 2/ ' 2 
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J 


47m’ e 


2 / 

/m; 


Kq and are zero- and first- 


order Bessel functions of the imaginary argument, respectively; and G and D are 
constants of integration. Four constants of integration (A, C, D, and G) are used, 
and four boundary conditions are required in order to obtain the dispersion equation 
that gives cu in terms of k. The continuity of the normal components of both u 
and E gives 


K’ ll( Kr o) 

K K i( K ’ r o) 


(ID 
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and 


c n’ J l( kr Q> 
n K^kr^ 


( 12 ) 


The continuity of the tangential component of E gives 


C = -A 


r^nmcuK 


K 




l( kr o) 


K i( K ' r () I o( & o) + r I i( Kr o) K o( K ' r o) 


(13) 


The pressure balance equation for electrons requires some explanation. The 
equation of motion for electrons is 


nm 3^ = -neE - vP = n v^ 1 - vp 
dt *e e *e 


(14) 


where 4? is the electric potential defined by E = -y4>. Integration across a boundary 
gives rise to the continuity of -ne4> + p , which can be written in the conventional 
notation 



(15) 


In equation (15), n is the average electron density on the boundary surface and is 
represented by the arithmetic mean of n and n’ . This boundary condition and 
equations (11), (12), and (13) yield the dispersion equation 


(16) 


where I Q = I Q (Kr 0 ); 1 ± = I^Kr^; K Q = K 0 (K’r 0 ); and Kj = K^K’r^. 

The right-hand side of equation (16) is a positive quantity; that is, co is real 
for all wave numbers k. Therefore, no instability occurs in the equilibrium con- 
figuration considered. 


2 

GO = 


zrne 2 (n + n’^kr^K^K’r^ 


K " m [ I o( kr Q) K l( kr Q) + v MS’SN)] [ I 0 K 1 * I : o K o 
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DISCUSSION 




The physical cause of stability can be understood by examining figures 2(a) to 
2(c). Figure 2(a) shows the equilibrium configuration, and figures 2(b) and 2(c) show 
the 1=1 and 1 = 0 perturbations, respectively. From figures 2(b) and 2(c), it is 
clear that on the concave side the electric lines of force M crowd" together and that on 
the convex side the lines of force are separated from each other. Thus, additional 
forces are produced that evidently decrease the bending illustrated in figures 2(b) and 
2(c). Mathematically stated, the magnitude of the electric field that results from the 
charge on the beam varies inversely with the distance of the electric field outside the 
beam. A perturbation that moves the beam surface inward provides greater electric 
force than a perturbation that moves the beam surface outward (fig. 2(b)). That is, 
since the direction of the electric force is outward from the beam, the resultant of the 
forces on the concave and convex sides in figure 2(b) is a force that decreases the 
bending. 


Electric lines of force 


t t t t > 

’ ' ' r Jr w w 


(a) Equilibrium. 



Electric lines of force 




(c) The 1=0 perturbation. 


Figure 2. - A column of the charged beam in equilibrium and perturbed. The 
arrows indicate the electric lines of force. 
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CONCLUSION 


From analysis of the stability of an electron beam injected into space, with the 
beam assumed to be confined to a cylindrical region and assumed to be in equilibrium 
under its electric charge and its thermal pressure, it is concluded that the beam 
charge produces a stabilizing effect. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January 27, 1970 
879-10-01-01-72 
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